home *** CD-ROM | disk | FTP | other *** search
/ Sprite 1984 - 1993 / Sprite 1984 - 1993.iso / src / lib / c / etc / symm.md / isinf.c.old < prev    next >
Text File  |  1990-08-10  |  2KB  |  100 lines
  1. /* 
  2.  * isinf.c --
  3.  *
  4.  *    Machine-dependent procedure to determine whether a double is 
  5.  *    infinity.
  6.  *
  7.  * Copyright 1989 Regents of the University of California
  8.  * Permission to use, copy, modify, and distribute this
  9.  * software and its documentation for any purpose and without
  10.  * fee is hereby granted, provided that the above copyright
  11.  * notice appear in all copies.  The University of California
  12.  * makes no representations about the suitability of this
  13.  * software for any purpose.  It is provided "as is" without
  14.  * express or implied warranty.
  15.  */
  16.  
  17. #ifndef lint
  18. static char rcsid[] = "$Header: /sprite/lib/forms/RCS/proto.c,v 1.2 89/01/07 04:12:18 rab Exp $ SPRITE (Berkeley)";
  19. #endif /* not lint */
  20.  
  21.  
  22. /*
  23.  *----------------------------------------------------------------------
  24.  *
  25.  * isnan --
  26.  *
  27.  *    Return whether a double is equivalent to infinity.
  28.  *
  29.  * Results:
  30.  *    1 if the number is infinity, 0 otherwise.
  31.  *
  32.  * Side effects:
  33.  *    None.
  34.  *
  35.  *----------------------------------------------------------------------
  36.  */
  37.  
  38. int
  39. isinf(value)
  40.     double value;
  41. {
  42.     union {
  43.     double d;
  44.     long l[2];
  45.     } u;
  46.  
  47.     /*
  48.      * Put the value into a union so we can check out the bits.
  49.      */
  50.     u.d = value;
  51.  
  52.  
  53.     /*
  54.      * An IEEE Std 754 double precision floating point number
  55.      * has the following format:
  56.      *
  57.      *      1  bit       -- sign of Mantissa
  58.      *      11 bits      -- exponent
  59.      *      52 bits      -- Mantissa
  60.      *
  61.      * If the exponent has all bits set, the value is not a 
  62.      * real number.
  63.      *
  64.      * If the Mantissa is zero then the value is infinity, which
  65.      * is the result of division by zero, or overflow.
  66.      *
  67.      * If the Mantissa is non-zero the value is not a number (NaN).
  68.      * NaN can be generated by dividing zero by itself, taking the
  69.      * logarithm of a negative number, etc.
  70.      */
  71.  
  72.     /*
  73.      * check the exponent
  74.      */
  75.     if ((u.l[0] & 0x7ff00000) == 0x7ff00000) {
  76.  
  77.     /*
  78.      * See if the Mantissa is zero.
  79.      */
  80.     if ((u.l[0] & ~0xfff00000) == 0 && u.l[1] == 0) {
  81.         /*
  82.          * Infinity.
  83.          */
  84.         return (1);
  85.     } else {
  86.         /*
  87.          * NaN.
  88.          */
  89.         return(0);
  90.     }
  91.     } else {
  92.     /*
  93.      * Normal.
  94.      */
  95.     return (0);
  96.     }
  97. }
  98.  
  99.  
  100.